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We write down the basic equations of Fermi-liquid theory for mixtures of fermions and bosons, 
an example being 3 He- 4 He mixtures at low temperatures. Basically the theory is identical to the 
one derived by Khalatnikov, but it is derived in a different way, and includes more discussion. A 
simplifying transformation of the equations is found where the coupling of the normal and superfiuid 
components appears in a simple form. The boundary conditions are discussed. 
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I. INTRODUCTION 

The Fermi liquid theory formulated by Landau has become a paradigm of what can be the effect of interactions 
in a Fermi systemi. Landau formulated the theory originally for liquid 3 He. Khalatnikov generalized this theory to 
mixtures of fermions and bosons, and applied it to mixtures of 3 He and 4 He @]. The purpose of this article is to 
reformulate Khalatnikov's theory. Although the theory is basically the same, our approach is different. We avoid 
several complications by concentrating on the laboratory frame, by directly starting with the linearized theory, and by 
using the osmotic energy. We construct a transformation that eliminates the coupling between the superfiuid velocity 
and the quasiparticle momentum. We discuss boundary conditions. We also try to interpret the theory by discussing 
ON ' how the momentum of a quasiparticle can be divided into three contributions. The purpose is to formulate general 
Fermi-Bose liquid equations that are needed to calculate the response on a vibrating wire 3 -^. 

A general introduction to Fermi liquid theory can be found in Refs. ^-(^. Although our emphasis is on subjects not 
found in these reviews, we try to be self-contained. Note that Khalatnikov's theory^ as well as the present theory 
make no assumption about the diluteness of the fermion component in the mixture, which is assumed in several other 
treatments***. Thus the Fermi-liquid theory of a pure fermion system can be obtained as a limiting case of the 
theory. 

We start in Sec. ILT1 with some basic definitions and introduce the noninteracting system. The basic assumptions of 
the Fermi-Bose liquid theory are given in Sec. IIIII In Sec. IIVI we define the phenomenological parameters that enter 
the theory. The equations of motion are formulated in Sec. [Vj The conserved currents are identified in Sec. I VI I and 
, some discussion is given. In Sec. IVIII we make a change to new variables where the coupling between the normal 
(— | ■ and superfiuid components appear only through their densities. The scattering of quasiparticles in the bulk and at 
surfaces is discussed in Sees. IVILT1 and HXl The hydrodynamic limit of the theory is discussed in Sec. IXl 

O ' 

II. BASIC DEFINITIONS 

> . 

. An important thing to realize is that the mass current density J is the same as the momentum density P. The 
equivalence J = P follows because the momentum of a particle is mass times the velocity, p — mv, i.e. the mass m is 
transported at velocity v. This relation is valid in condensed matter under standard (nonrelativistic) conditions since 
the changes of momentum and mass associated with the interaction field are negligible. This is valid also in crystalline 
material, where even in the idealized limit of infinitely rigid lattice, one should allow part of the momentum or mass 
current carried by the lattice. The equivalence J = P is essential in connection of equations (TlT)|) - (j2"Tj) below. 

We consider a system of particles having one type of fermions and one type of bosons, with masses tof and me. In 
this section, we describe the system in the absence of interactions. The free fermions have momenta p and energies 
e p = p 2 /2mp. The state of the fermion system is described by distribution function n p that takes values and 1 for 
each momentum state. The number density of fermions np is given by 



n F = J n p dr, (1) 

where dr = 2d 3 p/(2irh) 3 . The fermions have a spin but since we are not considering spin-depedent phenomena, it 
only appears as a factor of 2 in dr. The ground state consists of a Fermi sphere, which has filled states, n p = 1, for 
momenta p < pp and empty states, n p = 0, for momenta p > pp. The Fermi momentum pp is determined by the 
number density np — p f /3tt 2 H 3 . The momentum density of fermions is Jp = J pn p dr. The bosons are assumed to be 
condensed to a state with velocity i>b and all excited states of the boson system are neglected. The number density of 
bosons is denoted by Ub and the momentum density Jb = wb^bI'b • The total momentum density of both fermions 
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and bosons J = Jf + Jb is thus given by 

J = m B riBVB + J pn p d,T. (2) 



III. BASIC ASSUMPTIONS 



We now turn to the discussion of the interacting system. Landau's Fermi-liquid theory applies to low energy 
excitations of the system. The basic assumption is that when interactions are turned on, the low-energy part of the 
excitation spectrum remains qualitatively the same as it is in the non-interacting system. More precisely, one assumes 
that the momenta of the excitations remain the same when the interactions are turned on at constant densities of both 
fermions and bosons. The energies of the excitations can be shifted, but the equilibrium Fermi surface is assumed to 
remain unchanged. Also the quasiparticle energy is assumed to be linear in p close to the Fermi surface. Thus the 
excited states can still be specified by a quasiparticle distribution function n p , and the ground state corresponds to 
the filled Fermi sphere. 

A particular consequence of the assumptions is that Eq. (JTJ) remains valid in the interacting system. For the total 
momentum of an arbitrary interacting state we write, following Khalatnikov 2 

J = msnBVs + J pn p dr. (3) 

This is the same as Eq. ^ except that instead of the velocity of bosons it defines the super fluid velocity v s . 

In order to stress the non-triviality of Eq. (J3|) we mention that it is incorrect to deduce that the first term would 
be the momentum density of bosons and the latter that of fermions. The correct decomposition will be given later in 
equations ((22]) and ([24"]) . Note that an interacting system does not have a single boson velocity vb, and therefore it 
cannot appear in Eq. ([3]). 

In addition we use in the following the spherical symmetry of the system, the conservation laws of particle number, 
momentum and energy, and some estimates of orders of magnitude. 



IV. PARAMETERIZATION 

Above we used variables n pi ub and v s to specify the state of the system. In particular, the energy density could 
be written E({n p }, tib, v s ). Here the curly brackets indicate that E is a functional of n p , i.e. it depends on n p at all 
values of p. With respect to the variable ub, it is more convenient to change to the "osmotic energy" defined by 

E({n p },^ Bl v s ) = E({n p }, n B , v s ) - n B ^B - rap/4 . (4) 

Here the term — t£bMb effects the standard Legendre transformation from the density «b to the chemical potential 
fiB = dE({n p },nB,v s )/dnB- The term — riF/4 looks formally similar for the fermions, but is different because the 
chemical potential of fermions in a non-equilibrium state is not defined. Instead, we define a reference state. It has the 
equilibrium distribution n p ^ — 0(pf — p), where Q(x) is the step function. This equilibrium is taken to correspond 
chemical potentials j«g and /4° , superfluid velocity v s = and temperature T = 0. Thus the osmotic energy d3| 
still depends on the distribution function n p , and the effect of the — npfip^ term is merely a shift of energy so that 

quasiparticle energies are counted from /4°^ . 
The quasiparticle energy is defined as 



_ 5E({n p },fi B ,Vs) 



(5) 



p 



where the functional derivative is interpreted as SE = J e p 5n p dT. The linearity of the quasiparticle energy near the 
Fermi surface is satisfied by the choice 

e p = v F (p - p F ) + Se p . (6) 

The coefficient of p — pp defines the Fermi velocity v-p. Writing vp — pp/m* it defines the effective mass to*. The 
second term Se p in ^ can be written as 

5e p = (1 + a) + Dp ■ v s + J f(p,pi)(n pi - n < £ > )dr 1 , (7) 
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where SfiB — fJ-B — Mb ' ^ s ^ ne deviation of the boson chemical potential from its equilibrium value. We see that 8e p 
is linear in the deviation from equilibrium, and thus appears less important than the first term in ©. It was an 
important observation of Landau that this term still has to be kept in order to make a consistent theory. Being linear 
in the deviations of /ib > v s and n p from equilibrium, we see that Se p ([7|) has the most general form that is allowed by 
symmetry. For example the most general function of p that is linear in v s and does not depend on any other direction, 
has the form D(p)p ■ v s with some function D(p). For the same reason f(p,pi) cannot depend on the directions of p 
and p\ separately but only on their scalar product p ■ p\. 

We can further simplify 5e p ([7]). Since we are considering only excitations near the Fermi surface, we can approximate 
p w pfp in the Dp ■ v s term. Also, we estimate that f(p,Pi) changes essentially only on the momentum scale of pp. 
Since n p — rip^ is nonzero only in a thin shell around the Fermi surface, we can neglect the dependence of f(p,pi) 
on the magnitudes p and p\. Similarly, we can assume a and D as constants not depending on p. Thus Se p depends 
only on the momentum direction p, not on the magnitude p. We can now define an energy integrated distribution 
function 

,(o) 



<pp - n p ')vFdp, (8) 

where the argument of n p is p — pp. We get that / PlPl depends only on the scalar product p ■ p±. Thus we can 
expand f PtPl using Legendre polynomials Pi (x) 

* oo 

f(p ■ p') = ^f(P •p')=E FiPi (p ■ #)> ( 9 ) 



1=0 



where Pq(x) — 1, Pi{x) — x, etc and m*pY /ir 2 h 3 is the quasiparticle density of states at the Fermi surface. As a 
result we can write ([7]) into the form 

oo 

Se p = {l + a)S f iB + Dp F p ■ v s + ^ F t (P l (p ■ p')<j> p , ) p , , (10) 

1=0 

where (. . .) p denotes the average over the unit sphere of p. 

We have defined parameters to*, a, D and Fi's in equations ^ and (flU]) . There is one constraint between these 
required by translational invariance. The relation is conveniently stated by expressing D in terms of other parameters, 



to* \ 3 

This relation will be justified below in connection of equation ([23"]). Note that because of using the osmotic energy 
right from the start, we need to consider only one set of Fi parameters, in contrast to the three sets used in Ref. [2|. 
We also need to consider changes in the boson density. We parameterize 

5n B = Us — n-n' = -(1 + a)Sn F H ^<5/i B , (12) 

TObS 

where Sn-p — np — iip' . The equality of the coefficient 1 + a with the one in (JT]) and (IT0|) follows because they are 
second partial derivatives of i? (g]), 

d 2 E d^F driB 



dn Fd , B djUB 0n F = 1 + a - (13) 

In the limit of vanishing concentration of fermions, a reduces to the parameter a used in Refs. 8.9. In the same limit, 
the parameter s in (112[) reduces to the velocity of sound [as will be evident from equations (|18p and (|36[) below] . 



V. EQUATIONS OF MOTION 

In general the distribution function n p depends on location and time, n p (r,t), and similarly for other variables 
/ib (*",£) and v s (r,t). The quasiparticle distribution obeys the kinetic equation 
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where I p is the collision term. The fermion number and momentum conservation in collisions requires that 

J Ipdr = 0, J plpdr = 0. (15) 
Assuming small deviation from equilibrium, we linearize the kinetic equation (|14[) and get 

-£ +v F p-V [n p -^-5e p ) = I p , (16) 



where the derivative dn p ^ /de p should be evaluated at the unperturbed energy e p ^ = vp(p — pp). 
The superfluid velocity v s appearing in ([3]) is assumed to be curl freeii 

Vxb s = 0. (17) 

It is assumed to obey the ideal-fluid equation of motion, 

^ + — V MB =0. (18) 
at itib 

VI. CURRENTS 

Based on the equations above one should be able to derive conservation laws. In particular, one should be able to 
derive separate continuity equations for the two components 

to f ^ + V-J f =0, (19) 
at 

m B^ + V-J B = 0, (20) 
at 

and the momentum conservation 

dJ «■ 

- + v.n = o. (2i) 

where J = J F + Jb. In the following we derive linearized expressions for the fermion mass current J F , for the boson 
mass current Jb and for the momentum flux tensor II. 

We integrate the kinetic equation (jl~6|) over all momenta. Using ([lj and ((TO)) and comparing to ((1~9| we get 

J F = Dm F n F v s + ^(1 + -fi) [p(n p - n p 0) )dr. (22) 
to* 3 J 

From equations ([3]) and (|22j) we can eliminate the integral term and get 

( Dm*np\ to* , . 

J = m B n B - ! v s + J F . (23) 

V 1 + iFi/ TO F (l + ^l) 

This can be used to justify relation (TTTT) . Consider equilibrium but the whole system moving at constant velocity v. 
Because it has to be that J — (to f tt, f + 7n B n B )v : v s — v and Jp = TO F n F i», we get (jlip . 
Using © and ([2^)1 we calculate Jb = J J F and get 

Jb = (tobtt-b — Dmpnp)v s + D J p (n p — n p °^)dT. (24) 

Let us consider Eqs. (f2"2"j) and (IM1) in the frame where d s = 0. From Eq. (|2"4"]l we see that a single quasiparticle of 
momentum p has fraction Z) of the momentum carried by bosons. As the principal fermion forming the quasiparticle 
travels with velocity v = p/m*, it carries the fraction to f /to* of the quasiparticle momentum. The rest fraction of 
the momentum mpFi/Sm* is carried by other fermions. The total fermion contribution is visible in the second term 
of J F 
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In order to clarify the nature of the quasiparticle, we consider a simple model. Assume that the principal fermion 
forms a sphere of radius a and the bosons and other fermions are modeled as an incompressible ideal fluid of density 
p. In this model one can calculate m* = mp + 2ira 3 p/3, i.e. the effective mass is the fermion mass plus the fluid mass 
corresponding to half of the sphere. This result is obtained by calculating the momentum p = m*v from Newton's 
second law when a force is applied to accelerate the sphere, or, equivalently, by calculating the kinetic energy E/- 
and writing Ek — ^m*v 2 . This model gives a simple picture how the effective mass is increased due to the medium 
dragged along by the principal fermion. This is sometimes called "backflow" but this may be misleading since the 
momentum is in the "forward" direction (as m* > mp). 

We note that when considering a single excited quasiparticle in a finite container, one must also allow a compensating 
flow. Using the simple quasiparticle model as discussed above, the total momentum of an incompressible fluid in finite 
stationary container must vanish. Thus the momentum p = m*v of a quasiparticle must be compensated by an 
opposite flow of the medium. This compensating momentum arises from force applied by the walls of the container, 
and thus is separate from the momentum of the quasiparticle. 

For further insight, consider the case that the fermions are in equilibrium but moving with velocity v, i.e. Jp = 
mpnpv. We assume v s = 0. The mass density dragged in such a normal current is called normal fluid density p ni 
i.e. J = p n v. The complementary fraction of the total density p = mpnp + ™b"b that remains at rest, is called 
superfluid density p s = p — p n . From (|23[) we get 

Dm*np m*n-p 
Ps = m B n B - Y^T~f\/3' Pn = Y+F\Ji' (25) 

These expressions have a straightforward interpretation in terms of the momentum division introduced above. Namely 
Jp = mp(l + ^Fi)J/m* and J B = DJ . Solving J from the former and using Jp — mpnpv one gets the normal fluid 
density (|25|) . One also sees that the mass density Dm*np/{1 + -Fi/3) subtracted from the boson density in p s ((25)) 
just corresponds to the boson fraction that is bound to the quasiparticles. 

Next we check the momentum conservation (|2"Tj) . We take the time derivative of J ([3]), and use the kinetic equation 
of quasiparticles (|16p and the equation of motion of the Bose component (fTS|) . Using fj 15[) one finds the momentum 
conservation with the momentum flux tensor 

S = P(°> 1 + n B Sp B 1 + -L [pp (n p - <> - ^Se p ) dr. (26) 



p 

Here the constant term arising from the equilibrium pressure P^ ' is added. 

VII. ELIMINATION OF SUPERFLUID- VELOCITY COUPLING 

Our purpose in this section is to rewrite the theory so that the coupling between the superfluid and normal fluid 
components appears in a minimal way. For that purpose we define new quantities Sn p and c5e p , and write the kinetic 
equation in the form 

Q ( _ d n (0) 



. )f .5n p + ^ 5e P ) + v fP ■ V5n p = I p , (27) 



so that Se p will only depend on 6p B and on 5n p but not on v s . In order to achieve this, we scalar multiply the 
superfluid equation (|T5|) by Avppdn^ /de p and add it to the kinetic equation ([Tij]) . A is a constant that is fixed 
below. Comparing this to (1271) we find 

8n p = n p - — 5e p H Sp B ) , (28) 



Sir, = Sep H 5p B — Avpp ■ v s . (29) 

m B 



Similar to <\> p (8]), it is convenient to define 



%jj p = I 8n pp vpdp. (30) 
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where we have defined 



The next task is to express Sep in terms of 1pp. Doing this one notices that the dependence on v s drops out by 
choosing A — Dm* / (1 + i*i/3) and we get 

K °° Ft 

5i p = i Z f ^ B + 7^ — \—^( p i(P-P')^P')p'^ ( 31 ) 

1 + F ° ^ 1 + m+i F i 

m*D to* mp , . 

X = ; i r + 1 + a = = + 1 + a. 32) 

m B (l + ±Fi) m B (l + |Fi) to b V ; 

It is useful to express the densities and currents in terms of Sn p (|2"5)) or tpp ([3H)) . Here we give the fermion density 
Sii F = n F - «i 0) = 1 ( / <5n p dT - m ^q- ^Mb ) = oiJTTT f ^ _ ^Mb) , (33) 



1 + -Fb\./ tt 2 ^ 3 / 7T 2 ^ 3 (1 + F ) 

the fermion momentum (|22j) 

Jf = "~7 / ' P 5fi P dT = ^f^f (j3^p)p, (34) 

to* 7 7T^n. J 

and the momentum flux tensor (|26|) 

8 = P (0) t+ — J/z B Y+ — f ppSfipdr = P (0) Y+ — <5/i B Y + 3n F (ppi/;„)p. (35) 

TO B TO* 7 TO B 

In the momentum flux Q35p the second and third terms arise from the supcrfluid and the normal components, respec- 
tively. We also rewrite the boson conservation ([2"U| using (fT2")) as 

n B <9<fyt B ^ 9n F n mm 

^r + ^ v ^- mBi ^ =0 - (36) 

The equations (f2"7]) . and {jTSJ) constitute equations of motion in variables Sn p , <5/z B , and v s . 

We see that in the system of equations (|18D . (|2T|) . (|3~Tj) . and the coupling of the superfluid and normal 

components takes place through their densities via the terms proportional to K. Let us estimate this effect considering 
a problem of length scale a, frequency w, and relatively long mean free path I > a. We first consider fermion flow with 
velocity scale u, for which we estimate Snp ~ npu/vp. Assuming cua/s <C 1, the first term in (|36l) can be neglected, 
and it is important to balance the last two terms. We get that a fermion flow with velocity u generates a superfluid 
velocity of order K{np /np,)(Loa/vp)u. The opposite effect of v s ~ u on the quasiparticles can be estimated from the 
fermion conservation (TTSJ) and applying (j3"31 and We get that for oja/vp < 1 the induced fermion velocity is 

of order K{mB,/mp)(tjoa/vp) 2 u. For small frequencies wa/vp <C 1 the velocities induced in the other component are 
small, and the superfluid and quasiparticle components can be treated independently of each other. 

VIII. COLLISION TERM 

For detailed discussion of the collision term I p we refer to Sykes and Brooker— . In the collision of two quasiparticles, 
their total energy is conserved. It follows that the collision integral can be written as a function Sn p l|28p. The collision 
term can in a few cases be analyzed exactly^. A simpler approach is to use relaxation-time approximation. In its 
simplest version one assumes^ 



1 



C°)(eC°>+*ep-c-&.p)l, (37) 



where t is the relaxation time and coefficients c and b are chosen so that conditions (|15l) are satisfied. The relaxation 
time can also be specified by the mean free path I = vpr. The relaxation-time approximation leads to essential 
simplification because instead of using Sn p or Sn p , one can construct equations for %p p (|30p . which does not depend 
on the magnitude of the momentum. In the relaxation-time approximation the kinetic equation (1271) takes the form 

d 1 

— (tpp - Sep) + v F p ■ Vipp = -- (tpp - (tpp')p, ~ 3(p-pVp')p') • (38) 

It is also possible to introduce different relaxation times 77 (with I — 2,3,...) corresponding to spherical harmonic 
decomposition of tpp (|5"4l . 
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IX. BOUNDARY CONDITIONS 



The boundary conditions at walls depend on the scattering properties of the wall. Two common models for surface 
scattering are specular and diffuse scattering. We generalize the analysis of Ref. [l3| to Fermi-Bose liquids and include 
general interactions ([9]). We first consider the boundary conditions in the rest frame of a wall, and after that generalize 
to the case of a moving wall. 

We assume a wall with normal h that is inpenetrable to both bosons and fermions. From equation (|23[) we see 
that the vanishing of the particle fluxes J and J F normal to the wall imply vanishing of the normal component of the 
superfluid velocity, 

h-v s =0, (39) 
Specular scattering means that the quasiparticle momentum changes from p to p R = p — 2n(n • p), and thus 

n p = 7Jp_2fc(n-p)- (40) 

Before accepting this one should check that it conserves the quasiparticle energy ©, i.e. we calculate 

6e p - Se PR = Dp F (p - p R ) ■ v s + (F(p ■ p')<j) p , - F(p R ■ p')<f>p<) p > 
= 2Dp F {h ■ p)(n ■ v s ) + (F(p ■ p')<j>p - F(p R ■ P r )4> Pr ) p > 

= 2Dp F (n-p)(h -v s ) + (F(p ■ -4>p> R ])p> =0 (41) 

where we have used (fTUj) . p R ■ p' R = p ■ p', (|3"9")) and <j>p — <j)p' R — 0, which follows from P0|) . 

Diffuse scattering means that the quasiparticles reflected from the wall are in equilibrium evaluated at the exact 
quasiparticle energy ([6]) . Thus the distribution of the reflected quasiparticles is 

n p = n<°> (e^ + Se p - c) = + ^— {Se p - c) for n • p > 0, (42) 

where the constant c is determined by the condition that the particle flux (l3"4"|) to the wall vanishes, n ■ Jp — 0. More 
conveniently, the boundary is expressed as 

dn (0) _ 

8rip = c for n ■ p > (43) 

de p 

with a new constant c. The constant can be determined by evaluating the condition h ■ J F = for reflected quasipar- 
ticles, which allows to write the diffuse boundary condition as 



^ n p ~ IZ 1 — / h p' 5n p > dr' for n ■ p > 0. (44) 



4ir 2 h 3 dn^ 
m*p F de p Jn-p'KO 

The boundary conditions in the case of a moving wall can be found by changing the reference frame. In the 
laboratory frame the quasiparticle momentum p = p' + m F u, where p' is the momentum in the rest frame of the 
wall, and u is the velocity of the wall. (This relation follows because the momenta are assumed to remain unchanged 
when the interactions are switched on.) The quasiparticle distributions in the two frames are the same, n p — n' p> . A 
change arises in 5n p , 

X (0) / (0) X l dni0) - ( A r\ 

drip = n p -ny = n p , - n p / +mFU = Sn pl - -j—m F v F p ■ u (45) 

Based on this we find (\> p — (f> p + m F v F p ■ it, and using v s = v' s + u we get de p = Se'p + (to* — m F )v F p ■ u. The 
transformation for Sn p is then 

Sn p = Sn ' , — p F p ■ u (46) 

p de p 

which implies ip p — -0 p + PfP • u. 

Applying the transformation rules to (|3^|) . P0|) and pi| gives the superfluid condition 



h ■ v s — n ■ m, 



(47) 
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the specular boundary condition 



5n p = (5n p _2n(ftp) + — ; — 2p F n phu, (48) 



P 

and the diffusive boundary condition 

dn^de p f ,, dn(°) ,„ 2„. 

drip = 4— — - / n ■ p on p ' dr — p F \p + -i) u for n ■ p > 0. (49) 

dtp dT Jn p'<o d£p o 

In terms of ipg the specular and diffuse conditions can be written 

= ipp-2fi(n.p) + 2p F (n ■ p)(n ■ u). (50) 

2 

V>p out = -2(n • Pin^p in )p in + PF(Pout + g«) • m (51) 

with an average over half of the unit sphere (n • p; n < 0) . 

A consistency check for the boundary conditions is that the fermion particle current (f3~4"|) behaves as expected 

n ■ J F = m F n F n ■ u. (52) 

From the specular boundary condition (|50[) and from (|35p we deduce the vanishing of the transverse momentum flux, 

n-Sxn = 0. (53) 

X. HYDRODYNAMIC LIMIT 

The kinetic theory reduces to hydrodynamic theory in the limit of small mean free path of quasiparticles. We sketch 
the derivation here starting from the relaxation-time approximation. The quasiparticle distribution function ipp can 
be expanded in spherical harmonics 

oo / 

^ = EE c ^ m (p)- (54) 

1=0 m—-l 

The lowest order coefficients Cg and c™ can be related to fermion number (f3"3"|) and current (f3~4"|). Due to rapid 
scattering, all higher order coefficients are assumed to be small. We define normal fluid velocity v n by writing the 
current J F = n F v n . We multiply the kinetic equation (|3"5)l by pp — 1/3 and average over p. We get 

o -J 

-{(pp - 1/3) - Si*]) + M(PP - 1/3)P ' VVp) - -~{(PP - l/3)Vp)- (55) 

In the limit of small gradients and time derivatives, the time-derivative term and the coupling to third order term of 
([54| in the gradient term of f|55[) can be neglected. This allows to write the momentum flux tensor (|35|) into the form 



% = fl^ + ( ^Mb + 5P*) Sij - v ^ + -£--V. v nSij ) , (56) 



where effective pressure of the normal component 

1 



and the coefficient of viscosity 



SP* = KufSub + ^m*t|(l + F )5n F (57) 



V = -n F v F p F r. (58) 
5 
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The momentum conservation (|2ip combined with the flux (|56l) gives an equation of motion. Combined with the 
superfluid equation of motion (|18[) . the fermion conservation law (I19|) and the boson conservation law (|36[) gives a 
complete set of linearized hydrodynamic equations. Using variables (Snp, 6/j.b, v n ,v s ) the set can be written as 

^+n F V-v n = 0, (59) 

^^+PsV-v s+mB Kn F V-v n = 0, (60) 
s z at 

Pn ~dt +V5P * = V^Vn + ^VV -v n , (61) 

m B ^ + V^B - 0. (62) 

Here (I61[) is the Navier-Stokes equation describing the viscous normal component. 

One application of the equations (j59]) - (|62j) is to determine the velocities of first and second sound modes. The 
results are identical to those found by Khalatnikov 2 . 



XI. DISCUSSION 



For stationary phenomena, the Fermi-Bose liquid theory does not differ from the noninteracting system, when 
written in terms of proper variables. This follows because dropping the time derivative term in the kinetic equation 
(12711 . the only difference to noninteracting case is that the distribution 5n p is replaced by Sn p . The boundary conditions 
(P5 ]) - P9"j) and the observables (|3"3")) - (l3~4")) also are functions of Sn p and they depend on the interaction parameters to*, 
a, D and i*j's in a simple scaling manner, if at all. 

The force applied to slowly moving objects in a Fermi liquid in the ballistic limit was calculated in Ref. [3 All the 
results presented in Ref.lTil [Equations (13)-(19)] concern the time-indepent case. Therefore, their independence of the 
interactions parameters to*, a, D and i 7 } follows most simply from the general argument of the preceding paragraph. 

In time dependent problems, the interactions have important effect through the Se-tevm in the kinetic equation (|27l) . 
For small frequencies, however, the coupling to the superfluid motion can be neglected, as argued in Sec. I VIII Thus 
at low frequencies the response of Fermi-Bose liquid to external perturbation is the sum of independent superfluid 
and Fermi-liquid responses. 

The theory formulated here is applied to calculate the force on a vibrating cylinder in Refs. H, 0, Q, andHs! 
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